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As an extension of the former study on 2-dimensional systems, we simulate phase behavior of polymer-grafted col- 
loidal particles in 3 dimensions by molecular Monte Carlo technique in the canonical ensemble. We use a spherically 
symmetric square-step repulsive interaction potential, which has been obtained using self-consistent field calculation. 
In previous articles, we have studied these model colloids in 2 dimensions and found that these particles, although their 
interaction is purely repulsive, self-assemble into a string-like assembly, in a narrow region of physical parameter sets. 
In the present work, we show the existence of the string-like assembly in 3-dimensional systems and study the statistical 
properties of the arrangement of these strings. The average string length diverges around the region where the melting 
transition line and the percolation transition line cross, which has also been found in 2 dimensions. 

PACS numbers: 64.60.ah, 61.43.Er, 82.70.Dd 



I. INTRODUCTION 

We studied, in our previous works, the phase behavior of a 
polymer-grafted colloidal system in 2 dimensions 17 . In these 
studies, linear diblock copolymers were grafted onto the col- 
loidal particles, so that the hard sphere interaction between the 
particles are modified. 

We have derived a spherically symmetrical interparticle po- 
tential between a pair of the polymer-grafted colloidal parti- 
cles, using numerical self-consistent field (SCF) calculation. 
Our SCF results have shown that the pair interaction potential 
approaches, with some parameter sets, to repulsive step po- 
tential containing no attraction 1 . This potential was approxi- 
mated by square-step potential with a rigid core, 

<p(r) = oo r < cr\, 

<p{r) = e (> 0) <T\ < r < cr 2 , 
(j)( r ) = cr 2 < r, 

where r denotes the distance between the centers of the pair of 
the particles, cr\ denotes the diameter of the colloids, and <x 2 
is the diameter of the outer core formed by the polymer brush. 
A positive constant eo represents the repulsion originated from 
the grafted polymers. 

Utilizing molecular Monte Carlo simulation of the particles 
interacting via 0(r), we studied the phase behavior of the col- 
loidal system in 2 dimensions. At low temperature and high 
density, the particles, which have no attractive interaction, 
start to self-assemble and finally align in strings. We have 
called this effect "frogspawn effect" on account of the char- 
acteristic particle configuration. The string-like assembly is 
observed in small regions of the temperature and the density. 
The width of the potential step, 0-2/0-1, where the string-like 
assembly is found is limited to an interval 1.7 ^ 0-2/0-1 g 2.8, 
i.e. in the vicinity of 0-2/0-1 = 2.0. At this value 0-2/0-1 = 2.0, 
when the particles are closely packed on a straight line, both 
the hard inner cores, diameter o-\ , and the soft outer cores, di- 
ameter cr 2 , of the particles simultaneously pinch and confine 
each other. This means that the system tends to keep the lin- 
ear, i.e. the string-like, assembly and that the region of 0-2/0-1 



of the string-like assembly ranges around 0-2/0-1 = 2.0. 

We also discovered that the string-like assembly indicates 
a similarity to percolation transition and critical phenomena. 
The average string length diverges around the region where 
the melting transition line and the percolation transition line 
cross. The Fisher exponent, a critical exponent for this perco- 
lation transition, r = 1.9 is kept at any step width 0-2/ o-\. 

In the canonical ensemble (A^Vr-ensemble), the globally- 
isotropic string-like assembly was found as a metastable 
structure, i.e. an amorphous solid 1-9 , whereas the globally- 
anisotropic defect-free string-like assembly was discovered 
via 3-reservoirs method as the equilibrated structure 2 ' 3 . 

Malescio and Pellicane also simulated the same system at 
0-2/0-] = 2.0 and 2.5 and found dimers, trimers, and other 
various structures 8 ' 9 . The 2-dimensional system has recently 
been studied in experiments, where the string-like assembly 
has actually been confirmed 10 . These various assemblies were 
also found in the system at other step widths 11 , cr 2 /cri, and 
in a system composed of particles interacting via continuous 
repulsive potential 12 similar to </>(r). 

On the other hand, in 3 dimensions, a glass transition has 
been observed in the same model system 13 . A large number 
of ground states of the same model system have been dis- 
covered at zero temperature via genetic algorithms both in 2 
and 3 dimensions 1415 . However, the string-like assembly in 
3 dimensions was beyond scope of these recent works. De- 
spite a variety of findings in 2 dimensions shown above, the 
3-dimensional systems have been studied far less than the 2- 
dimensional systems. 

In the present work, we simulate our model system and 
study the phase behavior at finite T in 3 dimensions 1 . Sim- 
ulation methods are given in section II. Simulation results at 
C2/C1 = 2.0 are discussed in section III. At this step width, 
the string-like assembly tends to be stabilized 1 ' 4 . Therefore, 
at this cr 2 /cri, we check whether the string-like assembly also 
appears in 3 dimensions. Simulation results at other step 
widths are presented in section IV. 
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II. SIMULATION METHODS 

Monte Carlo simulation in the canonical ensemble (NVT- 
ensemble) via the standard Metropolis algorithm 16 ' 17 is per- 
formed in 3 dimensions, where N, V, and T denote the number 
of particles, the system volume, and the temperature respec- 
tively. Mersenne Twister is chosen as a random number gen- 
erator for our simulation 18 20 . The particles are, in the initial 
state, arranged on homogeneous fee lattices in a cubic sys- 
tem box with periodic boundary conditions. In one simulation 
step, a particle is picked at random and given a uniform ran- 
dom trial displacement within a cube of §Ao-\ long in each 
direction. A Monte Carlo step (MCS) is defined as N trial 
moves, during which each particle is chosen for the trial dis- 
placement once on average. After 1.0 x 10 6 MCS, by which 
the system relaxes to the equilibrium state except at low tem- 
perature and high density, we acquire data every 10 4 MCS and 
get 100 independent samples of particle configurations. The 
number of particles is fixed atN = 4000. 

0"i and eo are taken as the unit length and the unit en- 
ergy respectively. Dimensionless temperature is defined as 
k B T/eo, where k B T denotes thermal energy. The close-packed 
volume of the outer cores of the particles is denoted by Vo 
(Vo = No- 2 3 V2/2 in 3 dimensions). Dimensionless volume is 
defined as V/Vq. 

III. SIMULATION RESULTS AT <r 2 j(T X = 2.0 

Here we discuss simulation results at 02/0" 1 = 2.0, 
for which the string-like assembly was observed in the 2- 
dimensional case. 



A. Average string length at (t 2 I(T\ = 2.0 

The average string length at cr 2 /<Ti = 2.0 in 3 dimensions 
is given in Fig. 1 . The average string length has a sharp peak 
at V/Vq = 0.76 and k B T/eQ = 0.18, which means a divergence 
of the string length. Outside the vicinity of this peak, the long 
strings are absent. The position of the peak corresponds to 
the peak in the 2-dimensional system 1 , i.e. V/Vq = 0.64 and 
k B T/e Q = 0.12. 

Snapshots of the system taken at the peak are shown in 
Fig. 2. These snapshots indicate that the long strings form 
hexagonally-arranged cylindrical structures, which is similar 
to cylinder phase typically observed in diblock-copolymers. 

These results presented in this section demonstrate that long 
strings and the string-like assembly are also observed in 3 di- 
mensions, at 0"2/c"i = 2.0. 

B. Percolation phenomena at cr 2 /o"i = 2.0 

The occurrence probability of a percolated cluster at each 
parameter set, k B T/eo and V/Vq, is given in Fig. 3. A sharp 
boundary, i.e. percolation transition line, separates percola- 
tion and non-percolation regions. Compared with the results 



in 2 dimensions 14 , the percolated phase extends to high V/Vq, 
i.e. low density, regions. Since many degrees of freedom in 
high dimensions provide a large number of paths from one 
side of the system box to the opposite side, percolation thresh- 
old is lower in 3 dimensions than in 2 dimensions. 

Here we analyze, along the percolation transition line in 
Fig. 3, the cluster size distribution, n(s), where s denotes the 
cluster size. An example of n(s) on the percolation transi- 
tion line is presented in Fig. 4. This graph shows a relation, 
n{s) oc <r T , with r = 2.2. This t means the Fisher expo- 
nent, a critical exponent, of the percolation transition of our 
system. This power law and the Fisher exponent r = 2.2 are 
kept along the percolation transition line. This value, r = 2.2, 
in 3 dimensions is slightly higher than r = 1.9 in 2 dimen- 
sions 14 . Divergence of the string length in Fig. 1, similar to 
critical phenomena, is found around the percolation transition 
line. This result is the same as the 2-dimensional system. 



C. Thermodynamic phase diagram at 0-2/0; = 2.0 

Here we construct a thermodynamic phase diagram at 
o- 2 /o-i = 2.0 in 3 dimensions. 

In a region of V/Vq > 1.0 and extremely low temperature, 
our system is regarded as a system composed of hard parti- 
cles with diamger cr 2 , i.e. the outer cores 1 ' 4 . Crystalization of 
hard particle systems, which is called Alder transition, occurs 
at V/Vq ~ 1.5 in 3 dimensions 21 . Calculating the probabil- 
ity density of the square displacement of the particles in low 
temperature regions k B T '/e < 0.1 at V/Vq = 1.4 and 1.5, we 
confirm this Alder transition of the outer cores of our colloidal 
particles. The results, e.g. presented in Fig. 5, show a signif- 
icant change in this probability density between V/Vq =1.4 
and 1.5. Clearly, this change demonstrates the Alder transi- 
tion. 

In other regions of parameter space for V/Vq and k B T/eQ, 
we construct a thermodynamic phase diagram, utilizing the 
mean-square displacement of the particles (MSD) during 99 x 
10 4 MCS. Due to the Alder transition of the outer cores, MSD 
at V/Vq = 1.5 and k B T/£Q = 0.01 is set as the criterion to 
determine whether the system is in a fluid phase or in a solid 
phase. In the regions of the parameter space where MSD is 
smaller than this criterion, the system is considered a solid. 
On the other hand, when MSD is larger than the criterion, the 
system is a fluid or a coexisting phase of both the phases. 

The constructed thermodynamic phase diagram is shown 
in Fig. 6. At high k B T/eQ Alder transition of the inner cores 
of the colloids occurs, whereas Alder transition of the outer 
cores is not observed. At low k B T/e Q only Alder transition 
of outer cores is observed. At intermediate k B T/eo, the solid 
melts when not only V/Vq but also k B T /eo increases. The re- 
gion of the string-like assembly and the percolation transition, 
discussed in section III B, is located in the middle of the solid 
phase far from the two Alder transitions, i.e. crystalization. In 
addition to the simulation results in 2 dimensions, the present 
results in 3 dimensions show that the average string length di- 
verges around the region where the melting transition line and 
the percolation transition line cross 14 . 
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FIG. 1. The average string length at 0-2/cri = 2.0. (b) is the same graph as (a), taken around a peak of (a) with a fine resolution. Blue regions 
represent short strings and red ones long strings. Grids denote the points where the simulation is performed. Due to statistical errors, small 
humps are found in a region of V/Vo = 0.5 and k B T/eo > 0.8. 




FIG. 2. Snapshots of the system at a 2 l<T { = 2.0, V/V Q = 0.76, k B T/e = 0.18, and 1.99 x 10 6 MCS. Yellow spheres represent the centers of 
the particles and blue lines denote networks of overlaps between the particles, (a): taken in the lateral direction of the strings, (b): top view of 
the strings. 



These results are qualitatively consistent with the results in IV. SIMULATION RESULTS AT VARIOUS STEP WIDTHS 

2 dimensions 1 ' 4 . 

We have confirmed the existence of the long strings and the 
string-like assembly in 3 dimensions at 1x2/0"! = 2.0. Here, 
simulating the system at various step widths, we determine a 
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FIG. 3. Occurrence probability of percolated clusters for the same system as Fig. 1. cri/cn = 2.0. (b) is the same graph as (a), taken with a 
fine resolution. Blue regions represent non-percolated phase and red ones percolated phase. 




FIG. 4. Cluster size distribution, n(s) = M(s)/N, at o^/cn = 2.0, 
V/Vo = 1.38, and k B T/€o = 0.5. s denotes the cluster size and M{s) 
is the number of clusters with the size s found in the system. This 
parameter set is located on the percolation transition line in Fig. 3. 
log(.s)-log(n(i)) graph is shown. A solid black line shows a linear 
fitting of these data by y = —r * x + b with t = 2.2 and b = -1.1. 



range of o^/cn in 3 dimensions, where the string-like assem- 
bly is observed. 

We have discovered 1,4 that the string-like assembly is ob- 
served in the vicinity of the step width o^/o-j = 2.0. At 
o"2/o"t = 2.0, when the particles are closely and linearly 
packed on a straight line, both the inner cores and the outer 
cores of the particles simultaneously pinch and confine each 
other 4 . This means that, in this particle configuration, the in- 
ner core is pinched by the nearest neighbors and that the outer 
core by the next nearest neighbors. Therefore, the string-like 
assembly is stable in the direction of the straight line. When 
these straight strings are closely arranged in the same direc- 
tion, the string-like assembly is also stable in the lateral direc- 
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FIG. 5. The probability density of the square diaplacement of the 
particles, W SD , during 1 x 10 5 MCS at o- 2 /o-] = 2.0 and k B T/e = 
0.01. The result at V/Vq = 1.4 is plotted by a green solid line and 
VI Vo = 1.5 by a black broken line. The number of analyzed particle 
configurations, i.e. the number of samples, is equal to 9. 



tion since an additional energy cost is required to move and 
force one particle of a string into the adjacent string. Thus, the 
system tends to keep the string-like assembly at (Tilery = 2.0. 
This results in the region of the string-like assembly rang- 
ing around o^/cn = 2.0. Our previous simulation results 
for 2-dimensional systems 1 ' 4 have indicated such a range in 
1.7 % 0-2/ a-\ g 2.8. On the other hand, in 3 dimensions, 
this interval of 1x2/0"! decreases since the additional degrees 
of freedom for the 3 -dimensional systems reduce the effect 
of pinch and confinement of the particles arranged along the 
straight line 4 . We simulate the 3-dimensional system and con- 
firm this result. 

Percolation phenomena and the thermodynamic phase be- 
havior of the system at various step widths are also studied in 
the present section. 
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FIG. 6. Thermodynamic phase diagram at o~ 2 /cr l = 2.0. Dark grey 
regions represent solid phases and light grey fluid or coexistence of 
both the phases. Alder transition points of the inner and the outer 
cores are pointed by arrows above the panel. 



A. Average string length at various step widths 



V. CONCLUSIONS 

Using particle Monte Carlo simulation, we have studied the 
phase behavior of 3-dimensional systems composed of the 
particles interacting via the square-step repulsive potential, 
4>{r). The string-like assembly has been confirmed in 3 dimen- 
sions, in a solid phase. The long strings compose hexagonally- 
arranged cylindrical structures, which has a similarity with 
cylinder phase of diblock-copolymers. The string-like assem- 
bly is related to the percolation and the critical phenomena, as 
was also observed in 2-dimensional systems 1 ' 4 . 

At o"2/o"i = 2.0, when the particles are closely and linearly 
packed on a straight line, both the inner cores and the outer 
cores of the particles simultaneously pinch and confine each 
other. This determines the optimum value of the interval of the 
step width around cr 2 /cri = 2.0, where the string-like assem- 
bly is observed. We have found that the corresponding interval 
for 2-dimansional systems 1 - 4 is 1.7 g o^/cn g 2.8, which be- 
comes significantly narrow in the case of 3-dimensional sys- 
tems. This is originated from the additional degrees of free- 
dom in 3-dimensional space, which reduce the effect of con- 
finement of the particles arranged along the strings. 

The Fisher exponent, a critical exponent for the percolation 
transition, of the system is t = 2.2, in 3 dimensions, at any 
step width 0-2/0-1. This is slightly higher than r = 1.9 ob- 
served in 2 dimensions 1 ' 4 . 



The average string length at 1x2/0"! = 1.1, 1.5, 1.9, 2.5, 
and 3.0 are given in Fig. 7. Long strings and the diver- 
gence of the string length are not observed at these values 
of <j2lcr\. The same quantity as Fig. 7(d), but with high 
resolution, is shown in Fig. 8, which illustrates that this re- 
sult is independent of the resolution of the data. The interval 
of the step width where the string-like assembly is observed, 
1.7 g 0-2/cri ^ 2.8 in 2 dimensions 1 ' 4 , is significantly re- 
duced for the present 3-dimensional case. The long strings 
and the divergence of the string length are found only in the 
close vicinity of o^/cn = 2.0 in our simulation. This result 
is consistent with the discussion given above in the present 
section. 



B. Percolation phenomena and thermodynamic phase 
diagrams at various step widths 

The occurrence probability of percolated clusters and ther- 
modynamic phase diagrams at cr 2 /cri = 1.1, 1.5, 1.9, 2.5, 
and 3.0 are shown in Figs. 9 and 10. These results are qual- 
itatively consistent with the results at cr 2 /a"i = 2.0, given in 
section III. The Fisher exponent, t = 2.2, of the percolation 
transition of our system is also unchanged at these values of 
0-2/0-!. These results indicate that the discussion on the re- 
sults at 0-2/0"! = 2.0 also apply to the results at other values 

of 0"2/o"i. 
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occurrence probability of percolated clusters for cr^lcry = 2.0 is given in Fig. 3. 
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FIG. 10. Thermodynamic phase diagrams for 1X2/0" i = (a): 1.1, (b): 1.5, (c): 1.9, (d): 2.5, and (e): 3.0. The corresponding thermodynamic 
phase diagram for o- 2 lcr l = 2.0 is given in Fig. 6. Dark grey regions represent solid phases and light grey fluid or coexistence of both the 
phases. Alder transition points of the inner and the outer cores are pointed by arrows above the panel. 



